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ABSTRACT 

We show that element diffusion can produce large fluctuations in the initial helium 
abundance of stars. Diffusion time-scale, which in stellar cores is much larger than 
the Hubble time, can fall below 10 8 years in the neutral gas clouds of stellar mass, 
dominated by collisionless dark matter or with dynamically important radiation or 
magnetic pressure. Helium diffusion may therefore explain the recent observations of 
globular clusters, which are inconsistent with initially homogeneous helium distribu- 
tion. 
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1 INTRODUCTION 

The standard cosmological theory predicts that the epoch 
of nucleosynthesis ended with extremely homogeneous ele- 
ment distribution. Presently observed local abundance vari- 
ations are commonly attributed to stellar nuclear reactions 
and supernovae explosions. Yet such variations can be pro- 
duced by another process, namely element diffusion. The 
diffusion of elements in stars has been widely studied (e.g. 
Noerdlinger 1977; Proffitt & Michaud 1991; Chaboyer et al. 
1992; Richard et al. 2002). In stellar cores, where the diffu- 
sion time-scale is typically of order 10 11 years, diffusion can 
increase helium abundance by up to a few per cent, thereby 
accelerating their evolution. However, as we show in this 
paper, diffusion can make much stronger impact on helium 
abundance of the diffuse gas clouds of stellar mass. There- 
fore, stars produced by collapse of such clouds, will have high 
helium abundance from the moment of their formation. 

Element diffusion is caused by pressure gradients which 
impart different accelerations to elements with different 
atomic mass. For diffusion to be fast, low friction (colli- 
sion rate) between the particles has to be accompanied by 
large pressure gradient. However, to sustain pressure gra- 
dient for sufficiently long time requires an opposing force, 
such as gravity. In a system composed by gas only, increas- 
ing gravitational force is possible only by increasing den- 
sity, which would simultaneously increase friction. One way 
out to avoid this problem is provided by the collisionless 
dark matter (DM), which contributes to gravitational force 
without increasing friction. Efficient diffusion is also pos- 
sible in systems with dynamically important magnetic or 
radiation pressure. There strong gradients in the individual 
pressure components (magnetic, radiation, thermal) may ex- 
ist even when the total pressure is constant, thus allowing 
to avoid strong gravity /high density requirement. Finally, a 



more exotic alternative for efficient diffusion comes from the 
MOND model (Milgrom 1983; Bekenstein 2004), which pre- 
dicts stronger gravitational force for the low-mass rarefied 
clouds than the standard newtonian theory. 

We have found that prior to the end of reionization, dif- 
fusion should have produced many helium rich gas clouds of 
stellar mass that could make up to from several per cent (in 
the standard ACDM model) and to half of the total baryon 
budget (in the MOND model). Later, as the gas is accreted 
by small halos and recombines, diffusion may again become 
efficient on the stellar mass scales, until the ambient pressure 
and metallicity rise to sufficiently high values. Helium diffu- 
sion, therefore, may explain large fluctuations of primordial 
helium abundance inferred from the recent observations of 
globular clusters (Bedin et al. 2004; Norris 2004; Lee et al. 
2005; Piotto et al. 2005). 

In §2, 3 and 4 we study the diffusion of helium, respec- 
tively, in the DM dominated minihalos, in the MOND model 
and inside pressure- confined clouds. In §5 we consider the ef- 
fect of turbulence on the helium abundance variation. In §6 
we discuss the implications of our results. 



2 GRAVITATIONAL SEDIMENTATION 
INSIDE DARK MATTER MINIHALOS 

In absence of radiation and magnetic pressure, the dynamics 
of different atomic species is determined by gravitational 
field, partial pressure gradients and frictional forces due to 
collisions: 



g- 



V(mkT) ^V) - Vi 



(l) 



where Vi, a„ n-, and m-, are respectively the average veloc- 
ity, acceleration, number density and particle mass of the 
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species i, and is the friction coefficient between species i 
and j. Since Ty, which is in effect the time-scale for the par- 
ticle belonging to species i to be stopped by collisions with 
species j, is typically much smaller than all other relevant 
time-scales, the diffusion can be described as a steady state 
process with a\ — in the gas central mass frame. Further- 
more, when the abundances are still nearly homogeneous, 
we can take V{n\kT)/m = /xVp/p, where p and p are, re- 
spectively, the total gas pressure and density, and p is the 
mean molecular weight. Using this in eq. (f) we get 



j 



limn/mi 



Vp 



For maxwellian velocity distribution rij is given by 



3v27r /mi(mi +mj) 



mjfcT 



(2) 



(3) 



where <Tjj is the momentum transfer cross-section. For the 
scattering of neutral hydrogen and helium atoms we adopt 
the experimental value, cthoH = 0.9 x 10~ 15 cm (Khouw, 
Morgan & Schiff f968). For the ionized gas the cross-sections 
are typically much larger, so in this paper we restrict the 
calculations to the neutral medium. 

For hydrostatic equilibrium, Vp/p = g, and the stan- 
dard abundances (X — 0.75, Y = 0.25), the relative velocity 
of helium and hydrogen atoms is 



Va c - Vh = 80 m/s x 



10- 



(^LJ)- 1 (JL 

'■ J U0 3 m- 3 / VlOK 



-1/2 



(4) 



Assuming for simplicity the spherical symmetry we can de- 
rive the evolution of particles abundances from the continu- 
ity equation 



d t rii 



-r 2 V(r 2 niVi) 



(5) 



For systems with uniform temperatures combining eqs. 
(4) and (5) and replacing g with —GM/r 2 , where M is the 
total mass inside radius r, gives the characteristic time-scale 
for helium enrichment 



THc 



d t ln(n H e/n H ) 



-1I/2 



(6) 



where / g = p ga s/(pgas + Pdm) is the local gas mass fraction. 
Conveniently th c is independent of both mass and the radius 
of the system, but is a function of temperature and the gas 
fraction only. 

When THe is much larger than the age of the halo the 
helium abundance in the minihalos grows linearly with time 



AY 
~Y~ 



t 

THc 



= 6-10" 



3/ -l T -l/2 ^ 



1 + Z 

10 



-3/2 



t = teubble 
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(i 



3/2' 



0.15 



-1/2 



(7) 



where zt is the formation redshift of the halo and 

l_ + £\ 3 / 2 " 
+ Zi) 



/l+£\" 3 / 2 I" / l+z X 3 / 2 ' 

I 10 J Vl + Zf) 



= 5.4 ■ 10 s yrs x 
0.15 



-1/2 



(8) 



is its age at redshift z. Inside large halos whose virial tem- 



perature, Tvir, is far above the background gas tempera- 
ture, Tb g , the gas fraction is close to the mean cosmic value 
~ 0.15, so that AY/Y <JC 1. On the other hand in minihalos 
with T v i r < 2bg the gas accretion is suppressed by pressure. 
At virialization the dark matter density contrast , whose evo- 
lution is nearly independent of gas dynamics, reaches ~ 200, 
so that the gas fraction drops roughly to 0.15/200 ~ 10~ 3 . 
Moreover the gas expansion continues until Tb g falls below 
T v i r , thus further lowering / g . 

After the gas thermally decouples from radiation at z ~ 
150, its temperature drops adiabatically (Peebles 1993). 



T b g = 0.017(1 + z) 2 



(9) 



Later between z ~ 10 and z ~ 20 the gas is reheated by 
X-rays and the UV photons, followed by reionization. From 
eq. (6) we see that the optimal conditions for diffusion in 
DM minihalos are during the epoch when the temperature 
drops to its minimum, which is roughly between 2 and 10 
K. At these temperatures and for / g = 10 -3 , the change in 
helium abundance, AY/Y, is of order unity (eq. (7)). 

It should be noted, however, that AY/Y can not be 
large when T v i r <§C Tb g . The approximation we used so far 
for all elements, (im p g — X7(nikT)/m, will not be valid after 
diffusion had enough time to make AY/Y close to its equi- 
librium value (around (mEe/fBH- l)T v i r /pTbg = 2.5T v i r /Tb g 
for Tvir <S Tbg). So AY/Y can be of order unity only when 

T/ir ^ -*bg- 

The condition T v i r ~ Tb g , required for efficient diffusion, 
sets a characteristic scale for the individual gas cloud that 
can become helium rich. The virial temperature of a halo of 
mass M and density p is 



= 15 K 



M 



\ 10 3 M G 



2/3 



1 + Zf 
10 



(10) 



Setting Tvir ~ Tb g and using eq. (7), we find the mass of the 
gas cloud whose helium fraction is increased by AY/Y 



f s M = 0.2M e 



l + z 
10 



Hi 



l + z 



+ Zf 



3/2 



\l + zj 



3/2' 



0.15 



-1/2 



(11) 



Thus the increase of order unity in helium abundance can 
be achieved in the gas clouds of stellar mass. However, since 
the efficient diffusion is restricted to the regions with low gas 
fraction, prior to reionization no more than a few per cent of 
the baryons can reside in the helium rich clouds. In regions, 
whose evolution is still in the linear regime, the variation 
of helium abundance does not exceed 0.1-1 % (Medvigy & 
Loeb 2001). 

To verify that the assumption of hydrostatic equilibrium 
is a valid approximation, we need to check that the sound 
crossing time, t s , in a cloud is smaller than its dynamical 
time, tdyn = (Cp) -1 / 2 , while both of them are smaller than 
the Hubble time. The second condition, td yn < teubbie, is al- 
ways satisfied inside the virialized dark matter halos, since 
idyn/teubbie = (87r/3(5 d m) 1/2 ~ 0.2, where <5 dm ~ 200 is the 
dark matter density contrast. The ratio t s /td yn in the adia- 
batically compressed gas clouds is 



tdy 



= 0.6(1 + 5 e 



1-2/3 



0.2M 



1/3 



10 



-1/2 
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V 200 ) 



1/2 



0.02 



-1/3 



(12) 



where (5 gas is the gas density contrast. Thus inside the stellar 
mass clouds the hydrostatic equilibrium is indeed a reason- 
able approximation. 



3 DIFFUSION IN MOND MODEL 

Modified gravity has been suggested as an alternative to the 
collisionless dark matter (Milgrom 1983; Bekenstein 2004). 
In MOND theory in the limit of weak gravity, gravitational 
acceleration is given by g — ^fgWgo, where <;n = GM/r 2 
is the standard Newtonian value. Based on the galactic 
rotation curves the value of go is expected to be around 
10~ 10 m- s~ 2 . Using this in eqs. (4) and (5) we derive the 
new diffusion time-scale 

ll T ,l/2/£N sl/2 



W^T 1 '^^) 1 ' Vs = 0.8 • 10 s yrs x 
So 



T 1/2 (l + <5 gas ) 



1/3 



0.2M 

Mr, 



1/6 



0.02 



(13) 



The dynamical time for modified gravity is td yn 
(Gpg/gn) -1 ^ 2 , so the ratio t s /td y n in this case is 



ts , i/ 3 /0.2M 

_= 0.4(1 + * gas ) f_ 



1/4 



(14) 



showing that the assumption of hydrostatic equilibrium is 
again reasonable for stellar mass clouds. 

Unlike the case with newtonian gravity, in MOND th c 
is no longer scale independent. The cloud virial temperature 
in the MOND regime (g < g ) is T vir w 10 K(M/M ) 1/2 , 
so prior to reionization, when Tb g falls to a few kelvins, the 
diffusion is again most efficient in clouds of stellar mass. 
However, in MOND all gas may undergo efficient diffusion 
(i.e., thc ~ ^Hubble) and not just a small fraction that resides 
within dark matter minihalos with M ~ Mj oans . 



4 DIFFUSION IN PRESSURE-CONFINED 
CLOUDS 

Radiation pressure from the resonant Lya photons is capa- 
ble of supporting large gas clouds (Braun & Dekel 1989). 
Since hydrogen resonance photons are far more abundant 
than helium resonance photons, the radiation pressure ex- 
erts stronger force on hydrogen than helium atoms. The 
magnetic pressure presents a similar case. The magnetic field 
acts directly only on charged particles (i.e., electrons and 
ions). Those in turn collide with the neutral atoms, transfer- 
ring to them the acquired momentum. Since the collisional 
cross-section of ions with hydrogen atoms is several times 
larger than with helium atoms (Osterbrock 1961; Krstic & 
Schultz 1999), the magnetic pressure also exerts stronger 
force on hydrogen atoms. The eq. (2) can thus be expanded 
to account for all pressure components 



Vhc — Vfe _ V(a th p th + QradPrad + ObPb) 



THeH 



(15) 



where p t h, Prad and pb are, respectively, thermal, radiation 
and magnetic pressure, a t h = 1 — (imp/mec ~ 0.7, and 



Q ra d and ob are typically close to unity. When the total 
pressure, p — p t h + Prad + Pb, is constant, but radiation or 
magnetic pressure have strong variation on the scale R, eq. 
(15) becomes 



Vhc - Vh _ (1 - Q th )p _ 03p 

THoH pR pR ' 

The time-scale for helium enrichment is then 
R 



(16) 



THe 



M cl 

0.2M© 



V Hc - U H 

2/3 



= 1.7 • 10 yrs x 

1/3 



p/k 



100 K • cm- 



V r— v 



5/6 



(17) 



From the above equation it follows that when pressure is 
fixed, the diffusion is most efficient at the highest tempera- 
ture possible for the neutral gas (T ~ 10 4 K). In the inter- 
stellar medium of the Milky Way, where the metallicity and 
pressure are high (p/k ~ 2- 10 4 K • cm~ 3 ), it generally takes 
around a few million years for a cloud with T ~ 10 4 K ei- 
ther to cool down and collapse into a star, or to be destroyed 
through shockheating or ionization. Since for a given pres- 
sure, t He is of order 10 years, the expected helium variation 
due to diffusion in stars that currently form in the Milky 
Way, typically should not exceed one per cent. By contrast, 
in the first halos with T v \ r ~ 10 4 K, which form out of the 
metal-poor gas, and where the pressure is significantly lower 
(p/k ~ 10 2 — 10 3 K • cm -3 ), the typical life-time of a cloud 
may be comparable to th c . Therefore, significant variation 
of initial helium abundances may be possible for the early 
generations of stars. 



5 TURBULENT MIXING 

Small-scale turbulent motion tends to reverse the process of 
helium segregation by mixing the gas from regions with dif- 
ferent abundances. When helium abundance gradient, W, 
exceeds (Vkc — VkVKurb-Rturb, where 14 urb and R tU ib arc 
the characteristic velocity and the length scale of the tur- 
bulence, turbulent mixing prevails over helium sedimenta- 
tion. In case both operate at the same time, then in a 
cloud of radius R c \ helium abundance may vary by at most 
AY/Y = (Vhc - VH)i?ci/Kurbi?turb. When turbulence be- 
comes strong after helium sedimentation had some time to 
act unopposed, then over time the initial abundance gradi- 
ent is diluted by a factor ~ e ~ VturbRturbt/R ci . By contrast, 
large-scale turbulent motion, with R c \ <g R tU ib, has no effect 
on the small-scale abundance fluctuations. Furthermore, if 
the cloud is magnetized, the magnetic stress may be able to 
suppress the turbulent mixing altogether. 

In general there is no reason to expect strong turbulent 
motion in the low density gas clouds that were first to form 
out of cosmological density perturbation field (§2 and 3). 
Strong turbulence may arise later, as these clouds are ac- 
creted by larger objects and become pressure-confined, but 
its length scales, amplitude and duration are very uncer- 
tain. However, as has been pointed out by Bruston et al. 
(1981), while the impact of turbulence on the abundance 
variation in the interstellar clouds is hard to calculate, ob- 
servations suggest that it should be limited. In particular 
turbulence failed to destroy an order of magnitude variation 
in the abundances of different isotopes of carbon (Encrenaz, 
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Falgarone & Lucas 1975; Goldsmith & Langer 1978) and the 
separation between the different types of dust grains (Car- 
rasco, Strom & Strom 1973). Therefore it seems likely that 
the fluctuations of helium abundance may survive as well. 



6 DISCUSSION 

We have shown that diffusion can significantly increase the 
helium abundance of protostellar clouds. Similarly diffusion 
may also produce spatial variation of the deuterium (Brus- 
ton et al. 1981) and lithium abundances. The amplitude 
of the variation depends on many factors (magnetic field, 
outside pressure, turbulence etc.) and in some systems the 
effect of diffusion may still be negligible. However, as demon- 
strated by our calculations, the conservation of primordial 
abundances until the onset of stellar nucleosynthesis, cannot 
in general be taken for granted. 

Recent observations strongly suggest that stars in at 
least some globular clusters have been formed with enhanced 
helium abundance (Bedin et al. 2004; Norris 2004; Lee et al. 
2005; Piotto et al. 2005). It has been even suggested that the 
initial abundance of helium may be the missing "second pa- 
rameter" (Caloi & D'Antona 2005). Several scenarios involv- 
ing early pollution by helium-rich but metal-poor stars have 
been suggested as an explanation. However, Bckki & Norris 
(2005) have shown that such pollution is very unlikely, un- 
less the helium rich gas comes from the outside source and 
is kept in place by high external pressure. Therefore, diffu- 
sion provides a strong alternative as an explanation to high 
helium abundance. 
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